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ABSTRACT 

The  significance  of  year-round  relationships  between  synoptic  storm 
activity  and  sea-surface  temperature  (SST)  over  the  North  Pacific  Ocean 
for  the  period  January  1969  through  December  1978,  were  examined  using 
cross-correlation  analysis.  The  possible  existence  of  a  seasonal  depen- 
dence in  the  SST-wind  relationships  was  examined  using  sub-sample  record 
cross-correlation  analysis  and  model  tests  with  both  a  constant  and  a 

variable  mixed-layer  depth,  respectively.  Wind  forcing  data  was  repre- 

3 
sented  by  (U^),  friction  velocity  cubed,  and  wind  stress  curl  (CURLyx), 

as  computed  from:  (a)  High-pass  filtered  wind  components  only  (periods 

less  than  ten  days),  (b)  High-  and  Low-pass  filtered  wind  components 

which  includes  the  interaction  terms  and  (c)  the  unfiltered  (total)  wind 

components. 

3        3 
The  correlation  maps  of  SST  with  both  (U^)H  and  (U^)HL,  calculated 

from  the  ten-year  record,  showed  large  areas  of  significant  negative 

values  in  the  midlatitude  North  Pacific  Ocean  when  the  SST  was  lagged  by 

3 
any  amount  from  zero-  to  three-months.  Lag  correlations  involving  (U^) 

from  (c)  and  (CURL-ri)  were  considerably  less  significant.  Zero  lag  cor- 
relations between  3/3t  SST  and  the  various  wind  parameters  generally 
substantiated  the  above  results  except  for  revealing  a  large  and  signi- 
ficant positive  correlation  between  3/8t  SST  and  the  (CURL7T)  (from  (c)) 
off  the  \iest   Coast  of  North  America. 

The  sub-sample  record  examination  resulted  in  generally  low  signifi- 
cance levels  and  no  statistical  confidence  as  to  the  existence  of  a 
seasonal  variation  in  the  relationships.  Model  tests  involving  a  pre- 
scribed, seasonally  varying,  mixed-layer  depth  showed  that  there  was  no 
detectable  seasonal  variation  in  the  SST-wind  relationships. 
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I.   INTRODUCTION 

This  study  examined  the  relationships  between  synoptic  storm  activity 
and  sea-surface  temperature  (SST)  over  the  North  Pacific  Ocean  from  20°  N 
to  60°  N  and  110°  W  to  130°  E  for  the  period  January  1969  to  December 
1978.  The  hypothesis  tested  is  that  the  monthly  mean  sea  surface  tempera- 
ture responds  to  atmospheric  wind  forcing,  as  manifested  in  storms,  and 
that  the  amount  of  response  is  proportional  to  the  amount  of  wind  forcing. 
The  concept  for  examining  storminess  comes  from  the  strong  relationship 
that  exists  between  synoptic  storms  and  air-sea  exchange  (Simpson  (1969), 
Fissel,  Pond  and  Miyake  C1977),  Elsberry  and  Camp  (1978))  and  previous 
work  accomplished  in  this  area  by  Little  (1980). 

The  major  goal  of  the  study  was  to  obtain  a  quantitative  measure  of 
the  statistical  significance  of  the  year-round  relationships  between  SST 
and  wind  forcing  found  by  Little  C1980).  This  was  done  using  the  method 
of  Sciremammano  (1979)  and  the  results  are  shown  in  Section  III.  A 
second  major  goal  of  this  study  was  to  investigate  the  possible  existence 
of  seasonal  dependence  in  the  SST-wind  relationships.  This  was  done 
using  two  different  methods  and  these  results  are  shown  in  Section  IV. 
A  third  objective  of  this  study  was  to  examine  how  well  the  SST-surface 
wind  relationships  can  be  seen  using  conventional  data  on  storm  tracks 
and  storm  frequencies.  The  results  of  this  examination  are  shown  in 
Section  IV. 

The  approach  and  expectations  of  this  study  were  that:  1)  Applica- 
tion of  the  method  of  Sciremammano  (1979),  for  the  determination  and 
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presentation  of  statistical  significance  between  two  time  series  of  geo- 
physical data,  would  produce  a  quantitative  measure  of  the  significance 
of  the  year-round  relationships  between  SST  and  the  chosen  wind  forcing 
parameters,  friction  velocity  and  wind  stress  curl.  It  was  expected 
that  significant  negative  correlations  would  occur  when  the  SST  lagged 
the  wind  parameters,  in  support  of  the  hypothesis  that  anomalous  stormi- 
ness  leads  to  subsequent  ocean  surface  cooling.  2)  Examination  into  the 
seasonal  dependence  of  relationships  between  SST,  3/8t  SST,  and  the  wind 
parameters  would,  perhaps,  reveal  stronger  correlations  and  greater 
significance  during  the  summer  months  when  the  mixed  layer  is  shallowest 
and  weaker  correlations  and  less  significance  during  the  winter  months 
when  the  mixed  layer  is  deepest.  3)  Further  insight  into  the  dynamics 
of  the  relationships  between  SST,  9/9t  SST,  and  the  wind  parameters  could 
be  obtained  by  examination  of  a  particular  month  or  series  of  months  in 
the  ten-year  record  when  the  correlations  were  especially  significant. 
The  data  sets  used  in  this  study  were  obtained  from  Little  (1980). 
The  atmospheric  data  were  prepared  by  Fleet  Numerical  Oceanography  Center 

CFNOC)  from  a  special  ten-year  series  of  six-hourly  surface  wind  analyses, 

3 
Friction  velocity  cubed  (U^)  and  the  vertical  component  of  the  curl  of 

the  wind  stress  CCURL-^t)  were  calculated  from  the  wind  data.  These  para- 
meters  were  chosen  because  of  the  direct  relationship  between  (U^)  and 
mechanical  mixing  of  the  sea  by  the  overlying  atmosphere  and  the  well- 
established  relationship  between  the  wind  stress  curl  and  Ekman  conver- 
gence  and  divergence.  Six-hourly  values  of  both  (U^)  and  curl  were  com- 
puted using  a)  High-pass  filtered  wind  components  only  (periods  of  less 


12 


than  ten  days),  b)  a  combination  of  High-  and  Low-pass  filtered  wind 
components  which  includes  the  interaction  terms,  and  c)  the  unfiltered 
(total)  wind  components.  Monthly  means,  long  term  (10  years)  monthly 
means,  and  monthly  anomalies  of  each  of  the  six  wind  parameters  were 
calculated  for  comparison  with  monthly  SST  anomalies.  The  six  wind 
parameter  and  SST  anomaly  fields  were  truncated  to  reduce  "noise" 
using  Empirical  Orthogonal  Function  (EOF)  analysis.  Normalized  cross 
correlations  of  the  resulting  seven  data  sets  were  utilized  to  examine 
the  relationships  between  monthly  anomalies  of  SST  and  synoptic  storm 
activity,  as  represented  by  the  wind  parameters. 
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II.   DATA  SET  DESCRIPTION 

This  section  describes  the  origin  and  preliminary  processing  that 
was  done  to  the  data  sets  used  in  this  study. 

A.  SEA  SURFACE  TEMPERATURE 

A  ten-year  time  series  of  monthly  sea-surface  temperature  (SST)  data, 
covering  the  period  January  1969  through  December  1978,  had  been  pre- 
viously obtained  from  Scripps  Institution  of  Oceanography  Climate  Research 
Group  through  the  generosity  of  Jerome  Namias  and  Robert  Born  (Little 
C1980)).  The  data  are  geographically  located  at  166  grid  points  (equally 
spaced  at  5°  latitude-longitude  intervals)  in  the  North  Pacific  Ocean, 
covering  the  area  from  20°  N  to  60°  N  and  110°  W  to  130°  E  (see  Figure  1 
Little  (1980)).  Each  monthly  temperature  value  was  converted  from 
Fahrenheit  to  Centigrade  and  processed  to  compute  the  anomaly  by  finding 
the  difference  between  the  monthly  value  and  the  ten-year  average  for  the 
same  grid  point  and  month  of  year.  A  120-month  record  of  SST  anomalies 
at  each  of  the  156  grid  points  was  thus  constructed. 

B.  FRICTION  VELOCITY  AND  SURFACE  STRESS 

Monthly  anomalies  of  six  different  surface  wind  parameters  involving 

3 
the  cube  of  the  friction  velocity  (U^),  and  the  curl  of  the  surface 

stress,  CCURL^t),  had  been  previously  obtained  from  Haney,  Risch  and 

Heise  (1981)  by  Little  (1980).  All  six  wind  parameters  were  calculated 

from  six-hourly  surface  wind  analyses  over  the  North  Pacific  Ocean 
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prepared  by  Fleet  Numerical  Oceanography  Center  (FNOC)  for  the  same  ten 
year  period,  January  1969  through  December  1978.  The  six  wind  para- 
meters were  calculated  as  follows. 

Synoptic  storm  activity  having  time  scales  less  than  ten  days  was 

extracted  from  the  time  series  of  surface  wind  at  each  grid  point  by  the 

3 
use  of  a  symmetric  high  pass  Gaussian  filter.  The  wind  parameters  (U* 

and  CURL-^t)  were  then  calculated  from  a)  the  High-pass  filtered  wind  com- 
ponents only,  b)  a  combination  of  High-  and  Low-pass  filtered  wind  compo- 
nents, and  c)  the  total  unfiltered  wind  components,  respectively.  Monthly 
means,  long  term  (.10  years)  monthly  means  and  finally  monthly  anomalies 
of  each  wind  parameter  were  constructed  and  interpolated  from  the  FNOC 
Polar  Stereographic  Grid  to  the  156  point  SST  grid. 

Little  (1980)  accomplished  further  analysis.  All  seven  anomaly  data 
sets  (the  SST  and  the  six  wind  parameters)  were  detrended  and  standardized 
to  remove  linear  trends  in  the  anomalies.  Empirical  Orthogonal  Function 
(EOF)  analysis  was  then  performed  on  all  seven  data  sets  in  an  effort  to 
enhance  the  "signal"  by  reducing  the  "noise"  using  the  technique  of 
Preisendorfer  and  Barnett  (1977).  A  varying  number  of  EOF  modes  were 
retained  for  each  of  the  data  sets  resulting  in  the  final  seven  data  sets 
being  truncated,  or  filtered  versions  of  the  original  data  sets.  The 
amount  of  total  variance  retained  in  the  data  sets  was  of  the  order  of 
52  to  75  percent  (see  Little  (1980)). 

Possible  relationships  between  the  data  sets  are  examined  in  the 
following  sections.  Section  III  examines  the  relationships  that  exist 
year-round  between  SST  and  each  of  the  six  wind  parameters.  Based  upon 


15 


the  results  of  that  analysis.  Section  IV  narrows  the  field  of  interest 

3 
down  to  SST,  9/3t  SST  and  (U^)  calculated  from  the  combination  of  High- 

and  Low-pass  filtered  winds,  and  investigates  the  possibility  of  seasonal 
^    dependence  in  the  relationships.  Finally,  in  Section  IV,  the  sequence 
•    September,  October,  November  1977  is  examined  in  an  effort  to  gain 

further  insight  into  the  dynamics  of  the  relationships  between  SST  and 

(U*)HL. 


I 
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III.   YEAR  ROUND  RELATIONSHIPS  BETWEEN  WIND  AND  SST 

A.   CALCULATIONS  PERFORMED 

1 .  Raw  Correlations 

In  order  to  investigate  relationships  between  the  wind  and  SST, 
raw  correlation  values  were  determined  between  SST  and  each  of  the  six 
atmospheric  forcing  functions.  The  complete  data  records  of  120  months 
were  utilized  to  calculate  cross  correlations  between  SST  and  each  of 
the  wind  parameters  at  lags  of  minus  three-to  plus  five-months.  A  lag 
of  minus  three-months  means  the  wind  parameter  is  correlated  with  the 
SST  which  occurred  three  months  previous,  whereas  a  plus  three-month 
lag  means  the  wind  parameter  is  correlated  with  SST  occurring  three 
months  later.  Of  interest  in  this  study  is  atmospheric  forcing  of  the 
ocean  and  thus  emphasis  is  placed  on  evaluating  the  positive  lags. 

Correlation  values  were  computed  over  the  entire  grid  and  con- 
toured at  intervals  of  0.1.  Spatial  maps  of  the  raw  correlations  for 
minus  three-to  plus  three-month  lags  are  presented  in  Little  (1980) 
and  for  this  reason  have  been  omitted  here.  The  conclusions  by  Little 
(1980)  were  that  there  is  apparently  significant  negative  correlation 
between  the  SST  and  wind  forcing  based  on  the  (approximate)  signifi- 
cance level  criteria  used  in  his  study. 

2.  Significance  Levels 

Interpretation  of  the  raw  correlations  eventually  requires  some 
consideration  of  the  significance  levels.  Little  (1980)  made  use  of  a 
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critical  correlation  value  derived  from  Davis  (1975)  in  which  a  corre- 
lation, r,  was  considered  significantly  different  from  zero  if  |r|  > 
0.20.  The  critical  correlation  value  of  0.20  was  determined  by  a 
qualitative  estimate  of  the  number  of  independent  events  in  the  records. 
It  was  also  suggested  by  Little  (1980)  that  perhaps  a  quantitative  esti- 
mate of  the  critical  value  could  be  obtained  by  the  use  of  the  technique 
presented  by  Sciremammano  (1979).  Sciremammano  (1979)  derived  a  con- 
venient method  for  determining  and  presenting  significance  levels  of 
cross  correlations  which  is  equivalent  to  that  of  Davis  (1976).  This 
approach  normalizes  the  raw  correlations  by  a  large  lag  standard  error, 
a,  which  is  estimated  by  the  Bartlett  formula 

i  =+c» 

a^  =^  ^  CX(.i)  CY(i)  (1) 

i=_oo 

Here  CX(i)  and  CY(i)  are  the  discrete  autocorrelation  functions  of  the 

input  processes,  XCi)  =  X(iAt),  Y(i)  =  Y(iAt)  and  N  is  the  number  of 

discrete  observations. 

The  value  of  a   estimates  the  artificial  correlation  between  the 

two  parameters  which  arises  from  the  interplay  of  the  dominant  time 

scales  of  the  input  processes  (represented  by  the  summation  of  CX(i)CY(i) 

in  (J  II  and  the  record  length  N.  Normalization  of  the  raw  correlations 

by  a   immediately  shows  the  significance  of  the  correlations  since,  to 

a  good  approximation,  C99  =  2.6  a,  C95  =  2.0  a,  and  C90  =  1  .7  a  for  a 

time  series  involving  at  least  10  degrees  of  freedom.  As  shown  by 

Sciremammano  (1979),  this  requirement  is  fulfilled  when  n  ^  — t  _^  10. 

a 
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For  lagged  correlations  the  equation   is  modified  somewhat  and 
becomes 

i=+oo 


4  =  rn<    X^     CX(i)   CY(i)  (2) 


where  K  is  the  magnitude  of  the  lag  (K>0)  and  N-K  is  the  effective 
record  length.  Thus  the  only  difference  between  the  value  of  a   given 
by  (1)  and  (2)  is  due  to  the  slightly  shorter  record  that  is  necessary 
when  calculating  lagged  correlations.  This  dependence  on  the  effective 
record  length  is  not  very   strong  since  a^     is  inversely  proportional 
to  the  square  root  of  N-K,  and  usually  N>>K. 

Application  of  this  normalization  procedure  to  the  raw  correla- 
tions requires  a  practical  evaluation  of  the  summation  shown  in  (2). 
This  was  done  by  writing 


+1 

T 

'-1 


=  ^   CXCi)  CY(i)  (3) 


where  i^  is  a  bivariate  integral  time  scale  (expressed  nondimensionally 
in  units  of  At),  and  I  is  the  maximum  value  of  the  summation  index  i. 
Enfield  and  Allen  (1980)  suggest  the  summation  limit  I  be  determined 
such  that  successive  values  of  x,  obtained  by  increasing  I,  change  yery 
little.  For  this  study,  the  autocorrelation  functions  were  computed 
for  all  i  in  the  range  -23<i<23.  The  summation  limit,  I,  in  (3)  was 
determined  such  that  successive  values  of  x  changed  less  than  ^%.     This 
criteria  was  met  in  all  six  pairs  of  cross  correlations  (between  SST  and 


19 


each  wind  parameter)  for  I<14.  The  resulting  a   values  for  all  lags  were 
on  the  order  of  0.1  over  the  entire  grid. 

Figure  1  shows  a  plot  of  two  typical  autocorrelation  functions  at 
a  particular  gi^id  point  in  the  Mid-North  Pacific  Ocean.  It  is  seen  that 

the  time  scales  of  the  SST  events  are  on  the  order  of  eight  or  nine 

3 
months,  whereas  those  of  the  (.U^)HL  events  are  typically  two  to  four 

months.  Figure  2  shows  the  value  of  a   calculated  for  these  variables 
over  the  entire  grid. 

The  raw  lagged  correlations  (between  SST  and  a  given  wind  para- 
meter) at  each  grid  point  were  normalized  by  the  (appropriate)  a^  value 
for  that  grid  point,  yielding  normalized  lagged  correlations  on  the  order 
of  one  to  three.  The  results  are  presented  below.  As  noted  above,  a 
value  of  2  is  significant  at  the  95  percent  confidence  level  while  a 
value  of  3  is  significant  at  the  99  percent  level. 

B.  RESULTS 

3 
1 .  Lag  Correlations  between  (U^)  and  SST 

Figures  3  to  11  are  maps  of  normalized  correlations  between  the 

3 
three  (U^)  wind  parameters  and  SST  for  minus  three-to  plus  five-month 

lags.  Positive  lags  are  for  SST  lagging  the  wind.  Figures  3  to  5  contain 

3 
the  correlations  between  SST  and  (U^)H  computed  from  the  High-pass  wind 

components  only.  Figures  5  to  8  contain  the  correlations  between  SST  and 

3 
(II^)HL  computed  from  a  combination  of  High-  and  Low-pass  wind  components 

3 
and  Figures  9  to  11  contain  the  correlations  between  SST  and  (U^)TOT  com- 
puted from  the  unfiltered  (total)  wind  components.  In  each  map,  shaded 
areas  are  indicative  of  negative  correlations  and  the  contour  interval 
Is  1.0. 
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Figure  2.     Zero  lag  a  values  for  the  variables  SST  and   (U^)HL 

based  on  a  120-month  record.     Contour  interval    is  0.02. 
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Figure  3.  Normalized  correlation  between  (U^)H  and  SST  at  -3,  -2,  and 
-1  month  lag.  Contour  interval  is  1  and  several  closed 
contours  are  labeled.  Shaded  areas  denote  negative  values. 
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Figure  4.  Same  as  Fig.  3  except  for  0,  +1  and  +2  month  lags. 
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Figure  5.  Same  as  Fig.  3  except  for  +3,  +4  and  +5  month  lags. 
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Figare  6.     Same  as  Fig.  3  except  for  (U^)HL. 
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Figure  7.  Same  as  Fig.  6  except  for  0,  +1  and  +2  month  lags. 


27 


{U^)HL 


LAG  =  +3 


LAG  =  +4 


LAG  =  +5 


140E         160 


180 


160 


140         120W 


Figure  8.  Same  as  Fig.  5  except  for  +3,  +4  and  +5  month  lags 
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Figure  9.     Same  as  Fig.   3  except  for  (U*) TOT. 
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Figure  10.  Same  as  Fig.  9  except  for  0,  +1  and  +2  month  lags. 
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Figure  11.     Same  as  Fig.   9  except  for  +3,  +4  and  +5  month  1 


ags 
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Examination  of  the  map  series  reveals  negative  correlations 
over  large  areas  of  the  grid  for  the  zero  and  positive  lags.  This  result 
is  consistent  with  the  hypothesis  that  anomalous  surface  wind  forcing 
produces  subsequent  anomalous  surface  cooling.   The  significance  levels 
of  the  correlations  are  also  immediately  seen,  recalling  C95  =  2.0  a  and 
C99  =  2.6  a.  Large  areas  with  values  of  2.0  and  3.0  are  noted  in 
Figures  3  to  11;  however,  only  those  areas  which  are  clearly  greater 
than  S%   of  the  total  grid  area  will  be  considered. 

Figures  4  and  5  indicate  a  region  of  significant  negative  corre- 

3 
lation  in  the  Mid-North  Pacific  Ocean  between  (U^)H  and  SST  for  the 

zero-to  plus  three-month  lags.  The  area  is  located  at  approximately 

30°  N  180°  W  trending  northeast-southwest  and  consisting  of  normalized 

correlation  values  of  2.0  and  3.0.  A  similar  geographic  region  of 

significant  negative  normalized  correlations  is  seen  in  Figures  7  and  8 

3 
for  (U^,)HL  and  SST.  Here  also  the  region  appears  in  the  zero-to  plus 

three-month  lags.  Other  limited  areas  of  magnitude  2.0  and  3.0  are 

present  in  (nearly)  all  the  figures  but  are  not  considered  of  importance 

in  this  study  due  to  their  incoherent  spatial  distribution  and/or  the 

small  percentage  of  total  map  area  that  they  occupy.  The  presence  of 

3       3 
significant  correlations  of  (U^)H  and  (U^)HL  with  SST  indicates  that 

over  these  areas  of  the  North  Pacific  Ocean,  above  normal  synoptic  storm 
activity  during  a  month  occurs  with  and  is  followed  by  monthly  values 
of  below  normal  SST.  This  occurrence  at  lags  of  zero  out  to  plus  three- 
months  appears  to  show  that  the  ocean's  response  to  anomalous  storm 
activity  is  both  immediate  and  long  lasting. 
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A  brief  examination  of  the  minus  three-month  to  minus  one-month 

3        3 
lags  for  (U^)H  and  {U^)HL  (Figures  3  and  6)  indicates  an  area  of  signi- 
ficant positive  correlation  in  the  Gulf  of  Alaska  region.  This  appears 
to  show  ocean  forcing  atmosphere,  vice  atmosphere  forcing  ocean,  in  that 
warmer  than  normal  SST  is  correlated  with  subsequent  above  normal  stormi- 
ness  (and  cold  SST  with  below  normal  storminess).  This  result  will  not 
be  pursued  further  since  the  present  study  is  primarily  designed  (through 
the  particular  choice  of  atmospheric  parameters)  to  investigate  the 
possible  forcing  of  the  ocean  by  the  atmosphere. 
2.  Lag  Correlations  between  CURL-7T  and  SST 

Normalized  correlations  between  SST  and  the  three  wind  stress 
curl  parameters  for  lags  of  minus  three-to  plus  five-months  are  presented 

in  Figures  12  to  20.  The  format  of  the  figures  is  similar  to  that  for 

3 

CU^f) .  SST  is  again  lagged  relative  to  curl  with  a  lag  of  plus  one-month 

showing  the  correlation  between  the  curl  field  and  the  SST  field  which 
occurred  one  month  later.  The  contour  interval  is  1.0  and  negative 
correlation  areas  are  shaded.  Similar  results  of  large  areas  of  negative 
correlations  were  expected  between  the  curl  parameters  and  SST  for  zero 
and  positive  lags  since  anomalous  values  of  wind  stress  curl  indicate 
anomalous  Ekman  pumping. 

Examination  of  the  curl  map  series  (Figures  12  to  20)  immediately 
shows  the  existence  of  less  spatial  coherence  for  the  curl  than  was 
seen  for  the  friction  velocity  (Figures  3  to  11).  This  is  especially 
true  for  the  storm  related  curl  parameters.  In  addition,  although  negative 
correlation  areas  are  present  at  positive  lag,  they  are  less  dominate  and 
cover  less  total  grid  area  than  for  the  friction  velocity.  Areas  of 
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Figure  12.     Same  as  Fig.   3  except  for  (CURL^xlH. 
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Figure  13.     Sane  as  Fig.   12  except  for  0,  +1   and  +2  month  lags 
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Figure  14.  Sane  as  Fig.  12  except  for  +3,  +4  and  +5  month  lags 
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Figure  15.     Same  as  Fig.   3  except  for  CCURL^xjHL. 
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Figure  16.     Same  as  Fig.   15  except  for  0,  +1   and  +2  month  lags 
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Figure  17.     Same  as   Fig.   15  except  for  +3,  +4  and  +5  month  lags 
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Figure  18.  Same  as  Fig.  3  except  for  (CURL^t)!^. 
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Figure  19.     Same  as   Fig.   18  except  for  0,  +1,  and  +2  month  lags 
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Figure  20.  Same  as  Fig.  18  except  for  +3,  +4  and  +5  month  lags 
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magnitude  2  and  3  occur  at  (nearly)  all  lags  for  all  three  curl  parameters 
with  the  most  significant  ones  being  those  in  the  Gulf  of  Alaska  and 
westward  along  the  Aleutian  Chain  in  Figure  15  (0  and  +1  lag).  This 
pattern  is  consistent  with  storm  induced  Ekman  pumping  cooling  the  SST. 
Another  area  which  may,  perhaps,  be  significant  is  seen  at  the  plus  one- 
and  plus  two-month  lags  in  Figure  19.  A  positive  correlation  area, 
trending  northeast-southwest,  occurs  off  the  North  American  coast  and 
indicates  that  positive  SST  anomalies  occur  after  above  normal  wind 
stress  curl.  This  sequence  is  contrary  to  that  expected  from  Eknan 
pumping. 

Examination  of  the  above  year  round  lag  correlations  between  the 

six  wind  parameters  and  SST  appears  to  indicate  that  in  the  midlatitude 

3        3 
North  Pacific  Ocean  there  is  significant  correlation  of  (U^)H  and  (U^)HL 

with  SST  for  lags  zero-to  plus  three-months. 
3.  Correlations  Involving  9/3t  SST 

In  an  effort  to  further  determine  the  data  sets  which  provide 
the  most  significant  correlations,  the  time  rate  of  change  of  SST, 
9/9t  SST,  over  a  two  month  period  was  calculated  along  with  a  two  month 
average  value  for  each  of  the  three  friction  velocity  and  wind  stress 
curl  parameters.  The  resulting  data  sets  (each  119  months  long)  were 
used  to  calculate  raw  cross  correlations  (at  zero  lag)  which,  in  turn, 
were  normalized  by  the  large  lag  standard  error  calculated  from  the 
appropriate  fields  using  (1).  Figure  21  shows  the  normalized  correla- 
tions between  3/9t  SST  and  the  three  wind  stress  curl  parameters  and 
Figure  22  shows  the  normalized  correlations  between  9/3t  SST  and  the 
three  (U^)  parameters.  Examination  of  Figure  21  reveals  a  rather  large 
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Figure  21.  Normalized  correlations  at  zero  lag  between  3/8t  SST 
and  {CmLf)H,    (CURL2^)HL  and  (CURLz^)TOT.  Contour 
interval  is  1  and  shaded  areas  denote  negative 
correlations. 
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Figure  22.     Same  as   Fig.   21    except  for   (u2)H,    {U^)HL  and   {[}l)TOJ . 
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region  of  strong  positive  correlation  between  9/3t  SST  and  (CURL^t)!^ 
trending  northeast-southwest  off  the  North  American  coast.  This  area 
appears  to  indicate  a  region  where  positive  wind  stress  curl  correlates 
with  positive  8/9t  SST  and  is  clearly  opposite  the  results  expected 
from  Ekman  dynamics. 

Figure  22  shows  large  areas  of  significant  correlations  between 

3       3 
3/9t  SST  and  both  (U^)H  and  (U^,)HL,  as  was  seen  in  the  lag  correlations 

with  SST.  Comparison  of  the  top  and  middle  maps  indicates  improved 

correlation  values  in  the  Gulf  of  Alaska  region  due  to  the  presence  of 

the  existing  mean  flow.  The  two  storm  related  parameters  clearly  corre- 

3 
late  better  with  d/dt   SST  than  the  total  (U^,)  parameter.  The  better  of 

the  two  parameters  is  considered  to  be  the  friction  velocity  computed 

from  High-  and  Low-pass  filtered  components  since  it  includes  the 

synoptic  storm  activity  in  the  presence  of  the  existing  mean  flow,  and 

3 
in  practice  is  nearly  5  times  larger  than  (LI^)H  (Haney,  Risch  and  Heise, 

1981).  For  this  reason,  the  (U^)HL,  SST  and  9/at  SST  data  sets  are  used 

in  the  next  section  to  investigate  the  possibility  of  a  seasonal 

dependence  in  the  air-ocean  relations. 
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IV.  SEASONAL  RELATIONSHIPS  BETWEEN  (U^)HL  AND  SST 

3 
Examination  of  the  year-round  relationships  between  (U^)HL  and  SST 

in  the  previous  section  showed  significant  correlations  for  the  zero-to 
plus  three-month  lags.  It  was  further  desired  to  determine  if  any 
seasonal  variation  existed  in  these  correlations.  For  example,  one 
might  expect  the  SST  anomaly  to  be  more  highly  correlated  with  (anoma- 
lous) atmospheric  forcing  in  late  summer,  when  the  mixed  layer  is 
shallowest,  than  in  late  winter,  when  it  is  deepest. 

Two  approaches  to  investigate  the  possible  existence  of  seasonally 

dependent  relationships  were  considered.  In  the  first  approach  (Section 

3 
A  below),  the  correlations  between  monthly  anomalies  of  SST  and  (U^)HL 

were  calculated  for  each  of  the  12  starting  months  of  January  to 

December  and  for  lags  of  minus  three-to  plus  five-months.  The  use  of 

correlations  that  depend  upon  the  starting  month  provides  information 

on  the  seasonal  variation  of  the  wind-SST  relationships,  while  the  use 

of  lag  correlations  provides  information  on  the  time  scales  of  the 

events.  As  suggested  above,  one  might  expect  larger  correlations  to 

occur  when  the  starting  months  are  the  summer  months  (June,  July,  August), 

During  this  time  of  the  year,  the  climatological  ocean  mixed-layer 

depth  is  considered  to  be  the  shallowest.  Anomalous  storm  activity 

could,  perhaps,  more  easily  bring  cooler  water  from  below  the  mixed  layer 

and  cause  cooler  SST  values.  On  the  other  hand,  one  might  expect  less 

significant  correlations  when  the  starting  months  are  the  winter  months 
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(December,  January,  February).  The  climatological  mixed-layer  depth 
is  considered  to  be  deepest  at  this  time  of  the  year  and  more  intense 
or  more  persistence  storminess  would,  perhaps,  be  required  to  bring 
cooler  water  from  below  the  mixed  layer.  The  lagged  cross  correlations 
are  first  examined  over  the  entire  grid  for  each  starting  month  and 
lag,  and  then  they  are  evaluated  over  a  particular  geographical  area. 

The  second  approach  to  examine  for  possible  seasonal  relationships 
(Section  B)  tests  the  hypothesis  given  by 

3/3t  SST  =  -  MCU^)/(agh^)  (4) 

where  a  is  the  thermal  expansion  coefficient,  g  is  gravity,  h  is  the 
mixed-layer  depth  and  M  is  a  constant  "tuning  factor".  Such  a 
hypothesis  can  be  derived  from  simple  mixed  layer  theory  if  the  effects 
of  surface  heating  and  cooling  are  neglected  (see,  for  example,  Denman 
(1973),  equation  C18)  and  (25)).  Two  tests  of  this  equation  on  a  monthly 
time  scale  are  made  using  the  3/8t  SST  and  (U^)HL  anomaly  data  sets. 
In  the  first  case,  h  is  taken  to  be  a  constant,  while  in  the  second  case 
it  is  prescribed  to  vary  regularly  each  year  over  the  annual  cycle. 
The  key  question  examined  in  Section  B  below  is  whether  the  use  of 
such  a  seasonally  varying  mixed-layer  depth  significantly  influences 
the  normalized  correlations  between  the  parameters  on  the  left  and  right 
hand  sides  of  (4) . 
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A.     SUBSAMPLE   RECORD  CORRELATIONS  BETWEEN   (U^)HL  AND  SST 
1 .     Calculations  Performed 

To   investigate  the  possibility  of  a  seasonal   dependence  in  the 

3 
air-ocean  relations,  the  complete  120-month  records  of  (LI^)HL  and  SST 

anomalies  were  divided  into  sub-sample  records.     Twelve  separate  records 
were  constructed  -  one  for  each  of  the  months  January  through   December  - 
with  each  record  containing  the  ten  anomaly  values   (one  from  each  year) 
for  that  month.     For  example,  the  first  sub-sampled  record  of  SST  con- 
tained the  SST  monthly  anomaly  values   for  January  1969,   January  1970, 

etc  to  January  1978.     The  second  sub-sampled  record  contained  the 

3 
monthly  anomaly  values  for  all   ten   February's,   and  so  on.     The  (U^)HL 

records  were  treated  similarly.     The  sampling  interval    in  the  sub-sample 

records  was  12  months,  whereas  the  interval    for  the  complete  record  was 

one  month. 

a.  Raw  Correlations 

3 
Raw  cross  correlation  values  between  SST  and  (LI^,)HL  were  com- 
puted from  each  of  the  12  monthly  records  for  lags  of  minus  three-to 
plus  five-months.  For  example,  a  zero  lagged  correlation  for  the  start- 
ing montK  of  January  includes  the  cross  correlations  obtained  from  the 
January  1969  SST  and  January  1969  (u2)HL;  January  1970  SST  and  January 
1970  (uJ)HL;  etc  to  January  1978  SST  and  January  1978  (uJ)HL.  The 
number  of  realizations  for  this  January  zero  lag  correlation  is  ten. 
When  considering  January  at  plus  five-months  lag,  for  example,  the  cross 
correlations  are  obtained  from  the  January  1969  (U;^)HL  and  June  1969  SST; 
January  1970  (U^)HL  and  June  1970  SST;  etc  to  January  1978  (uJ)HL  and 
June  1978  SST  and  the  number  of  realizations  is  ten.   For  negative 
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correlations  the  number  of  realizations  is  nine.  Thus  there  are  typically 
nine  or  ten  realizations  in  the  seasonal  lagged  correlations  compared 
to  115  to  120  realizations  for  the  lagged  correlations  based  on  the 
complete  record. 

The  seasonal  lagged  raw  cross  correlation  values  ranged 
from  0.2  to  0.8  with  a  large  percentage  of  them  between  0.5  and  0.7. 
This  can  be  compared  with  the  complete  record  raw  cross  correlations 
in  which  the  range  was  from  0.1  to  0.3  with  a  large  number  near  0.2. 
The  increase  is  largely  due  to  increased  artificial  correlation  as  is 
shown  below. 

b.  Significance  Levels 

The  question  of  statistical  significance  must  once  again  be 
considered  and  in  order  to  evalute  the  significance  levels,  seasonal 
values  of  the  large  lag  standard  errors  were  calculated  for  individual 
lags  and  starting  months.  The  o^   values  are  wery   dependent  upon  record 
length.  Reduction  of  the  sample  size  from  typically  120  for  the  com- 
plete record  to  typically  ten  for  the  sub-sample  record  causes  an 
increase  in  o^^   from  a  value  of  approximately  0.1  to  values  ranging  from 
0.3  to  0.4.  The  expected  artificial  correlation,  a,  increased  by  a 
factor  of  (approximately)  /Tz  =   3.46  as  the  record  length  (or  number  of 
realizations)  was  decreased  by  a  factor  of  (approximately)  12.  Evalua- 
tion of  the  bivariate  integral  time  scale,  x,  given  by  (3),  becomes 
more  difficult  due  to  the  less  well-behaved  nature  of  the  autocorrela- 
tion functions  which  were  computed  for  all  i  in  the  range  -23<i<23 
(Figure  23).  The  criteria  of  successive  values  of  f  changing  less  than 
^7o   as  the  summation  limit  I  was  increased  by  1,  could  not  be  achieved, 
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even  for  the  maximum  summation  limit  1=23.  For  this  reason,  the  large 
lag  standard  error  in  the  present  case  of  the  sub-sampled  record  was 
estimated  by  a^  =  /I2  a^,  where  a^   is  given  by  (2)  and  calculated  from 

K  K  i\ 

the  complete  record.  The  Qj,  values  thus  computed  were  (approximately) 
0.35  and  were  spatially  coherent  over  the  entire  grid  area.  As  an 
example,  Figure  24  shows  a  values  for  the  zero  lag  case. 

The  raw  lagged  correlations  at  each  grid  point  for  each 

* 
of  the  12  starting  months  were  normalized  by  the  a  value  for  that 

grid  point,  yielding  normalized  lagged  correlations  on  the  order  of 

1  to  2  as  shown  below. 

2.  Results 

3 
a.  Lag  Correlations  between  (U^)HL  and  SST 

Figures  25  to  28  are  maps  of  normalized  correlations  between 

3 
monthly  anomalies  of  CU^,)HL  and  SST  at  plus  one-month  lag  and  starting 

months  of  January  through  December.  The  correlations  for  the  starting 

months  of  January,  February,  and  March  are  shown  in  Figure  25,  while 

the  other  starting  months  are  in  Figures  26  to  28.  In  each  map,  shaded 

areas  are  indicative  of  negative  correlations  (e.g.  above  normal  stormi- 

3 
ness,  CU^)HL,  followed  one  month  later  by  below  normal  SST)  and  the 

contour  interval  is  1. 

Examination  of  the  map  series  reveals  negative  correlations 

over  large  areas  of  the  grid.  This  result  agrees  with  the  results 

from  the  complete  year-round  record  and  is  consistent  with  the  hypothesis 

that  anomalous  storminess  produces  subsequent  cooling.  The  significance 

levels,  however,  are  not  as  large  due  to  the  increased  value  of  o^   over 
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Figure  24.  Zero  lag  a     values  for  the  variables  SST  and  (U^)HL, 
Contour  interval  is  0,08. 
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Figure  25.  Subsample  record  correlations  (nondimensional ized  by 
o^)   between  SST  and  (u|)HL  for  +1  month  lag  and 
starting  months  January  (top),  February  (middle), 
and  March  (bottom).  Contour  interval  is  1  and  shading 
Indicates  negative  correlations. 
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Figure  26.  Same  as  Fig.  25  except  for  April,  May  and  June. 
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Figure  27.  Same  as  Fig.  25  except  for  July,  August  and  September, 
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Figure  28.  Same  as  Fig.  25  except  for  October,  November  and 
December. 
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G^,.     It  was  felt  that  the  existence  of  a  seasonality  in  the  relation- 

3 
ships  between  (U^)HL  and  SST  would  be  evidenced  in  the  map  series  by  one 

starting  month,  or  perhaps  a  series  of  starting  months  having  clearly 
higher  (nondimensional )  correlations  than  the  rest.  Figures  25  to  28 
do  not  show  such  a  result  and  appear  to  show  less  areas  of  magnitude  2 
and  3  and  thus  less  significance  than  the  complete  year-round  record. 
Further  investigation  into  the  seasonality  was  done  by  considering  a 
specific  portion  of  the  grid  or  an  index  area. 

b.  Central  North  Pacific  Ocean  Index  Area 

An  open  ocean  index  area  consisting  of  18  grid  points  cover- 
ing the  area  30°  N  to  40°  N  and  170°E  to  165°  W  was  chosen  to  further 
consider  the  question  of  seasonality  since  this  region  contained  the 
most  significant  correlations  based  on  the  complete  year-round  record. 
An  average  normalized  correlation  over  the  index  area  was  computed  from 
the  map  series  in  the  previous  section  for  each  starting  month  from 
January  through  December  and  all  lags  from  minus  three-to  plus  five- 
months  (Table  I).  The  values  in  Table  I  for  the  zero  and  positive  lags 
are  nearly  all  negative.  The  magnitudes,  however,  are  nearly  all  less 
than  2.  This  is  less  significant  than  the  average  value  (at  one  month 
lag)   of  -2.84  obtained  for  the  index  area  using  the  complete  year- 
round  record.  A  graphical  presentation  of  the  seasonal  correlations  for 
zero-,  plus  one-,  and  plus  two-month  lags  (Figure  29)  appears  to  indicate 
that  the  most  significant  [negative)  correlation  between  (LI^)HL  and  SST, 
averaged  over  the  index  area,  occurs  during  the  winter  months  December, 
January,  and  February,  whereas  the  least  significant  relationships  occur 
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Figure  29.     Plot  of  normalized  seasonal   correlations  between  the 
variables  SST  and   (U^)HL  averaged  over  the  Gulf  of 
Alaska   Index  Area.     Lags  of  zero   (top),   plus  one  (middle) 
and  plus  two   (bottom)  months  are  shown. 
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in  the  summer  months  of  June,  July  and  August.  This  is  different  from 
the  expected  result  noted  above. 

Because  of  the  generally  low  significance  of  these  correla- 
tions and  their  variation  during  the  year  and  the  fact  that  the  index 
area  was  chosen  a  priori  on  the  basis  of  significant  year-round  correla- 
tions, it  is  felt  that  no  firm  conclusion  can  be  drawn  from  the  results 
of  this  section.  That  is,  it  cannot  be  stated  with  any  statistical 
confidence  whether  or  not  there  is  a  seasonal  variation  in  the  lag  cor- 
relations. In  order  for  this  method  of  sub-sampled  record  analysis  to 
be  fruitful,  a  much  longer  total  record  length  is  required. 

B.  MODEL  TESTS 

1 .  Correlations 

The  second  approach  to  investigating  seasonal  relationships  be- 

3 
tween  CU*)HL  and  SST  was  one  of  a  test  of  the  hypothesis  given  by  (4), 

with  particular  attention  to  the  effect  of  a  seasonally  varying  mixed- 
layer  depth.  The  time  rate  of  change  of  SST  over  a  two  month  period  and 
the  average  CU^,)HL  over  the  same  two  month  period  were  calculated  from 
the  complete  120-month  record.  The  resulting  119-month  records  were 
used  to  compute  correlations  in  which  two  cases  were  considered. 

In  the  first  case,  the  mixed-layer  depth,  h,  was  held  constant 
over  the  entire  record.  Raw  cross  correlations  for  each  of  the  166  grid 
poiints  were  calculated  between  the  two  data  sets  and  normalized  using 
the  large  lag  standard  error,  a,  for  the  119-month  record  given  by  (1). 
Typical  values  of  a  were  0.07  to  0.09  with  all  values  less  than  0.1. 
The  resulting  normalized  correlations  ranged  in  magnitude  from  1  to  3, 
as  shown  previously  in  Figure  22  (middle)  and  discussed  further  below. 
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In  the  second  case,  the  mixed-layer  depth,  h(t),  is  prescribed 
to  vary  according  to 

27T(t-t^) 

h(t)  =  H  +  AH  cos  (— =|y— ^)  (5) 

where  H  =  75m  represents  an  annual  average  mixed  layer  depth,  AH  =  50  m 
represents  the  amplitude  of  the  seasonal  variation  and  t  its  ohase. 
In  (5),  t  runs  from  1  to  120  months.  This  function  allows  the  mixed- 
layer  depth  to  oscillate  each  year  from  a  maximum  of  125  m  in  month  t 
to  a  minimum  of  25  m  in  month  t=+6   For  example  if  t  =1,  the  maximum 

0  0 

depth  occurs  fn  Janiiajry' anrl  ralntrajujn  in  July  wh-ije.  If  t  =  2,  the  maximum 

occurs  in  February  and  minimum  in  August.  Raw  cross  correlations 

3      2 
between  3/9t  SST  and  (U^)HL/hCt)  were  computed  and  then  normalized 

using  the  large-lag  standard  error  calculated  as  in  (1).  Correlations 

were  calculated  for  all  values  of  the  phase,  t  ,  from  1  to  12  allowing 

the  deepest  mixed-layer  depth  to  occur  successively  in  each  of  the  12 

months. 

2.  Results 

The  middle  map  in  Figure  22,  presented  in  Section  III,  shows 

3 
the  normalized  correlations  between  8/9t  SST  and  (U^)HL.  This  repre- 
sents a  test  of  (4)  and  (5)  with  AH  =  0,  i.e.  constant  h.  Two  areas  of 
significant  correlation  (magnitude  2  and  3)  occur  in  the  figure  and 
appear  to  support  the  hypothesis  that  positive  values  of  anomalous 
storminess  are  well  correlated  with  negative  values  of  8/9t  SST.  The 
nine  grid  points  in  the  Gulf  of  Alaska  region  of  Figure  22,  which  have 
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normalized  correlation  values  of  magnitude  3  or  greater,  were  considered 
to  investigate  the  influence  of  prescribing  a  seasonally  varying  mixed- 
layer  depth. 

Figure  30  shows  the  normalized  correlation  value,  averaged  over 
the  Gulf  of  Alaska  Index  Area,  and  presented  as  a  function  of  the  phase 
of  the  mixed-layer  depth.  Figure  31  shows  a  similar  plot  for  the 
Central  North  Pacific  Ocean  Index  Area  used  earlier.  Examination  of 
Figures  30  and  31  reveals  that,  in  general,  yery   little  increase  in 
significance  is  obtained  by  the  introduction  of  a  seasonally  varying 
mixed-layer  depth.  For  both  index  areas,  it  appears  that  the  chances 
of  reducing  the  correlation  is  greater  than  the  chances  of  increasing 
it.  In  addition,  the  phase  of  the  mixed-layer  depth  that  results  in  the 
best  correlation  is  different  for  the  two  index  areas.  In  the  Gulf  of 
Alaska  Region  (.Figure  30)  the  best  (most  negative)  normalized  correla- 
tions occur  when  the  mixed-layer  depth  is  prescribed  to  be  deepest  in 
March  (and  shallowest  in  September),  while  the  poorest  (least  negative) 
correlations  occur  in  the  case  when  it  is  deepest  in  August.  In  the 
Central  North  Pacific  Ocean  Region  (Figure  31)  the  best  normalized  cor- 
relations occur  in  the  case  when  the  mixed-layer  depth  is  deepest  in 
December,  whereas  the  poorest  occur  in  the  case  when  it  is  deepest  in 
May  or  June.  In  both  of  these  areas,  the  introduction  of  a  seasonally 
varying  h  (even  if  optimized  as  to  phase)  does  not  improve  the  correla- 
tions.  It  thus  appears  that  while  there  is  (considerable)  support 
CFig.  22)  for  the  hypothesis  given  by  (4)  (i.e.  3/3t  SST  and  (uJ)HL 
negatively  correlated)  no  significant  improvement  is  obtained  by  the 
introduction  of  a  seasonally  varying  mixed-layer  depth  to  modulate  the 

(U^)HL  term, 
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In  order  to  give  additional  observational  support  to  the  (sta- 
tistically significant)  year-round  correlations  found  between  synoptic 
storm  activity,  represented  by  (U^)HL,  and  the  development  of  monthly 
SST  anomalies,  represented  by  3/8t  SST,  the  synoptic  developments 
observed  during  the  Winter  of  1977  are  described  in  the  next  section. 
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V.  SYNOPTIC  RELATIONSHIPS 


In  an  effort  to  gain  further  insight  into  the  dynamics  of  the  rela- 

3 
tionships  between  9/8t  SST  and  (U*)HL,  three  questions  were  considered. 

3 
First,  are  the  best  correlations  between  8/3t  SST  and  (U^)HL  occurring 

3 
only  when  (U^)HL  is  positive  or  only  when  it  is  negative,  or  are  the 

correlations  the  same  in  both  cases?  Second,  are  there  one  or  two  partic- 

3 
ular  cases  or  months  in  the  ten-year  time  series  when  SST  and  (U^)HL 

are  extremely  well  correlated  and  if  so,  can  this  high  correlation  be 

seen  by  comparing  the  SST  anomaly  maps  with  conventional  storm  track 

data  for  the  same  time  period?  Third,  are  the  best  correlations  in  the 

North  Pacific  Ocean  occurring  along  or  displaced  from  the  mean  climato- 

logical  storm  track? 

3 
In  considering  the  first  question  of  positive  or  negative  (U^)HL 

values  leading  to  the  best  correlations,  the  119-month  data  records  for 

3  3 

8/9t  SST  and  CU*)HL  were  utilized,  with  (U^,)HL  being  the  key.  For  each 

of  the  166  grid  points  and  all  t,  t  =  1  to  119,  it  was  determined  if 

3 
the  (U^,)HL  anomaly  at  t  was  greater  than  zero,  or  not.  The  positive 

3 
CU^)HL  anomalies  at  time  t  were  correlated  with  the  corresponding  SST 

3 
value  (Figure  32).  For  positive  (.U^)HL,  large  areas  of  negative  correla- 
tion occur,  indicating  above  normal  storminess  is  correlated  with 

3 
decreasing  SST  (3/9t  SST  <  0).  Similarly,  for  negative  (U^,)HL,  large 

areas  of  negative  correlation  also  occur  indicating  that  below  normal 

storminess  is  correlated  with  increasing  SST  (3/9t  SST  >  0).  Figure  32 
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(U^)HL  Positive 


(U^)HL  Negative 


140E        160  180  160  140         120W 


Figure  32.  Correlations  between  (U^)HL  and  d/dt   SST^for  positive 

Ctop)  and  negative  (bottom)  values  of  (U^)HL  separately, 

Contour  interval  is  0.1  and  shaded  areas  denote  nega- 
tive correlations. 
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appears  to  indicate  little  difference  in  the  correlations  among  the  two 
distributions.  Both  above  normal  and  below  normal  values  of  (U^)HL 
correlate  equally  with  8/3t  SST,  with  neither  being  clearly  dominate. 
None  of  the  correlations  in  Figure  32  are  yery   significant  on  a  purely 
statistical  basis  because  the  large  lag  standard  error  (artificial  cor- 
relation) for  these  data  sets  are  approximately  /2  times  that  of  the 
total  record.  Thus,  raw  correlations  must  be  near  0.3  to  be  signifi- 
cant at  the  95%  confidence  level. 

3 
To  determine  if  there  was  a  particular  period  in  which  SST  and  (U^)HL 

were  extremely  well  correlated,  the  Gulf  of  Alaska  Index  Area  described 
earlier  was  chosen  for  examination.  The  two  time  series  for  t  =  1  to 
119  were  averaged  over  the  index  area  and  plotted  in  Figure  33.  This 
figure,  as  much  as  any  other,  shows  the  short  time  scales  (approximately 
1  to  2  months)  which  lead  to  the  estimate  of  1/100  for  the  number  of 
independent  events  in  the  record.  It  is  because  of  the  small  a   values 
that  a  correlation  of  0.3,  while  small,  is  highly  significant.  The 
sequence  September,  October,  November  1977  was  chosen  for  further  exami- 
nation because  of  the  apparently  good  correlation  at  that  time.  The 
actual  CU^)HL  and  SST  anomaly  values  for  each  of  the  months  September, 
October,  and  November  were  plotted  and  compared  with  the  principal 
tracks  of  centers  of  cyclones  at  sea  level,  as  given  in  Mariners  Weather 

Log  (1978)  and  the  observed  climatology  of  principal  low  tracks  (Klein 

3 
1956).  This  was  done  in  an  attempt  to  see  if  the  anomalies  in  (U^)HL 

could  be  related  to  the  apparent  "storm  track  anomalies".  The  results 

are  shown  in  Figures  34,  35  and  36. 
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(U^)HL 
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Figure  34.  SST  anomaly  values  (°C),  (U^)HL  anomaly  values  (units  10 

(cm/sec)-^)  and  cyclone  tracks  (Mariners  Weather  Log)  for  the 
month  of  September  1977.  Contour  interval  is  0.5  and  shaded 
areas  denote  negative  values.  Heavy  dashes  denote  climato- 
logical  tracks  (Klein  1956). 
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Figure  35.     Same  as  Fig.   34  except  for  October,   1977. 
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Figure  36.  Same  as  Fig.  34  except  for  November,  1977. 
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3 
The  September  and  October  (U^)HL  maps  (middle  frame  in  Figures  34 

and  35)  both  show  a  large  region  of  strong  positive  anomaly  values  in 

the  Central  Pacific  Ocean  with  negative  anomaly  regions  on  either  side 

to  the  East  and  West.  This  agrees  quite  well  with  the  general  storm 

track  patterns  as  depicted  in  Mariners  Weather  Log  and  shown  in  the 

bottom  frame  of  Figures  34  and  35.  In  addition,  the  September  storm 

tracks  show  some  isolated  (and  obviously  anomalous)  tracks  in  the 

Central  Pacific  Ocean  which  appear  to  be  associated  with  a  sharp  trough, 

flanked  by  ridging  (i.e.  no  storms)  on  either  side.  The  October  pattern, 

on  the  other  hand,  shows  a  generally  southward  displacement  of  the 

storm  track  stretching  West  to  East  across  most  of  the  Pacific.  The 

3 
association  between  November  (U^)HL  anomaly  values  and  the  storm  tracks 

(Figure  36)  seem  to  agree  quite  well  in  certain  areas.  For  example, 

the  above  normal  values  in  the  Gulf  of  Alaska  region  with  negative  or 

below  normal  values  to  the  South  and  also  in  the  Far  Western  Pacific 

Ocean  north  of  Japan,  seem  to  match  the  storm  tracks.  But  there  is 

also  a  region  of  distinct  disagreement  in  the  Central  Pacific  Ocean 

between  30°  -  40°  N  where  the  storm  tracks  appear  numerous  while 

3 
(U^)HL  is  below  normal. 

In  considering  the  SST  response  to  atmospheric  storminess,  the 

absence  of  storm  activity  in  the  eastern  part  of  the  Pacific  Ocean  and 

the  Gulf  of  Alaska  region  during  the  month  of  September  appears  to  have 

brought  additional  warming  during  October.  As  the  October  storminess 

effect  began  to  be  felt  in  the  Gulf,  however,  the  trend  appears  to  have 

reversed,  reducing  the  November  SST  values  to  near  normal. 
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Considering  the  qualitative  nature  of  any  estimate  of  cyclone  fre- 
quency, it  is  perhaps  surprising  to  see  as  much  agreement  between  the 
surface  wind  parameter,  (U^,)HL,  and  an  estimate  of  cyclone  frequency 
anomalies  using  data  from  the  Mariners  Weather  Log.  Given  a  sufficient 
amount  of  data  over  the  ocean,  it  is  felt  that  the  two  should  tend  to 
agree  quite  well.  However,  (U^)HL  should  be  better  related  to  SST 
anomalies  because  of  its  direct  physical  connections  via  the  process 

of  turbulent  vertical  mixing.  While  it  can  be  seen  that  at  certain 

3 
times  and  places  (.U*)HL  can  be  associated  with  a  qualitative  estimate 

of  cyclone  frequency  anomalies,  the  comparison  is  yery   subjective. 

Finally,  Figure  22  (Section  III)  is  examined  to  determine  where  the 

best  correlations  occur  in  relation  to  the  mean  cl imatological  storm 

track.  The  large  areas  of  magnitude  2  and  3  to  the  north  in  the  Gulf 

of  Alaska  region  and  to  the  south  in  the  Central  Pacific  region  appear 

to  indicate  that  the  best  correlations  occur  to  the  right  and  left  of 

the  mean  track,  with  a  weak  minimum  in-between.  Thus  the  anomalous 

storm  activity  caused  by  North-South  displacements  or  meanders  in  the 

stom  tracks  seem  to  show  up  the  best  in  the  correlations  with  d/dt   SST, 
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VI.  CONCLUSIONS 

The  determination  and  presentation  of  significance  levels  of  cross- 
correlations  using  the  normalization  method  of  Sciremammano  (1979)  worked 
very   well  in  this  study.  For  the  atmospheric  and  oceanic  variables  con- 
sidered in  this  study,  the  year-round  data  sets  contained  approximately 
100  independent  events  (out  of  the  120  realizations)  and  the  large  lag 
standard  errors  were  on  the  order  of  0.1.  Attempts  at  sub-sample  record 
analysis  reduced  the  realizations  to  approximately  10  (i.e.  down  by  a 
factor  of  12)  and  increased  the  a   values  to  approximately  0.3  to  0.4. 

The  application  of  the  normalization  method  to  the  year-round  data 

sets  revealed  that  of  the  six  wind  parameters,  the  correlation  of  SST  with 

3  3 

(IJ^)H  and  with  ('J*)HL  showed  large  areas  of  significant  negative  values 

in  the  midlatitude  North  Pacific  Ocean  when  the  SST  was  lagged  by  any 

amount  from  zero-to  three-months.  In  addition,  similar  large  areas 

3 
appeared  in  the  Czero-lag)  correlations  between  SST  and  either  (U^)H  or 

3 
CU^)HL.  These  results  appear  to  clearly  support  the  hypothesis  that 

anomalous  synoptic  storm  activity  on  the  time  scale  of  a  month  or  more 

is  coincident  with  and  is  followed  by  ocean  surface  cooling.  The  most 

representative  data  set  is  felt  to  be  (U*)HL  since  it  considers  the 

synoptic  storm  activity  in  the  presence  of  the  existing  mean  flow,  and 

3 
the  fact  it  is  nearly  an  order  of  magnitude  larger  than  (U^,)H.  Zero 

3 
lag  correlations  between  SST  and  the  CU^)  parameters  showed  two  areas  of 

significant  correlation  in  the  North  Pacific  Ocean  and  support  the  above 

conclusions. 
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The  (zero-lag)  correlation  between  3/3t  SST  and  (CURL^t)!©!  revealed 
a  large  region  of  significantly  positive  values  off  the  North  American 
West  Coast  which  may  indicate  an  important  region  where  positive  wind 
stress  curl  anomalies  are  associated  with  monthly  SST  warming.  This 
result  is  opposite  the  correlation  expected  from  Ekman  dynamics. 

Examination  into  the  possible  existence  of  a  seasonal  dependence  in 

3 
the  relationships  between  SST  and  (U^)HL  gave  results  which  supported 

the  over-all  hypothesis  through  large  areas  of  negative  correlations 
but  did  not  show  values  which  were  as  significant  as  those  from  the 
year-round  record.  The  sub-sample  record  gave  generally  low  signifi- 
cance levels  and  no  statistical  confidence  as  to  the  existence  of  a 
seasonal  variation.  The  introduction  of  a  seasonally  varying  mixed- 
layer  depth  did  not  improve  the  correlations  over  what  was  obtained 
with  the  constant  h.  Thus,  the  present  data  sets  tend  to  suggest  that 
there  is  little  or  no  seasonal  variation  in  the  wind-SST  relationships. 

The  investigation  into  the  dynamics  of  the  relationships  between 

3 
SST  and  CU*)KL  revealed  that  the  correlations  were  about  the  same 

3 
regardless  of  whether  CU^t)HL  was  above,  normal  or  below  normal,  with 

neither  being  clearly  dominate.  The  sequence  September,  October, 

3 
November,  1977  showed  good  correlations  between  (U^)HL  and  3/8t  SST  but 

3 
attempts  to  see  if  anomalies  in  the  maps  of  (U*)HL  could  be  related  to 

"apparent  storm  anomalies"  in  the  Mariners  Weather  Log  was  only  partially 

successful.  At  certain  times  and  places  the  (U^)HL  parameters  could  be 

associated  with  a  qualitative  estimate  of  cyclone  frequency  anomalies 

but  the  comparison  was  yery   subjective.  It  was  found,  however,  that  the 
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best  correlations  occur  to  the  right  and  left  of  the  mean  cl imatological 
storm  track,  with  a  minimum  in-between. 

Uhile  this  study  has  revealed  large  areas  of  statistically  signifi- 
cant correlations  between  monthly  anomalies  of  synoptic  storm  activity 
and  SST,  the  actual  raw  correlations  were  not  yery   large.  Typical  raw 
correlations  were  0.3  with  the  large-lag  standard  error,  a,  near  0.1. 
This  generally  low  value  of  the  raw  correlation  is  an  indication  of  the 
importance  of  other  factors,  both  atmospheric  and  oceanic,  which  influence 
the  SST  but  which  were  not  considered  in  this  study. 

Suggestions  for  future  research  deal  with  an  evaluation  of  the 

3 
accuracy  with  which  (_U^)HL  depicts  storm  activity  and  the  possible 

relationships  between  CCURL7t)T0T  and  9/3t  SST  resulting  in  the  region 
of  significant  positive  correlation  observed  off  the  North  American 
Coast.  While  it  is  felt  that  CU^)HL  is  a  useful  parameter  in  that  it 
gives  a  quantitative  measure  of  synoptic  storm  activity  which  has  not 
been  previously  available,  it  may  be  desirable  to  check  its  accuracy 
in  depicting  storm  activity  by  making  this  same  type  of  synoptic  com- 
parison in  a  more  data  rich  region  such  as  the  continental  United  States. 
A  possible  hypothesis  for  the  significant  positive  correlation  region 
off  the  North  American  Coast  suggests  a  relationship  between  positive 
wind  stress  curl  and  below  normal  low  level  (stratus)  cloudiness.  The 
resultant  increased  solar  radiation  could  possibly  lead  to  positive  SST 
anomalies  at  one-  or  two-months  lag  if  this  process  were  strong  enough 
to  dominate  the  cooling  due  to  Ekman  pumping.  Cloud  cover  data  for  the 
area  and  time  period  would  be  required  to  investigate  such  a  relationship. 
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